pdf(file="Asymm.pdf")
hist(ztrace[ztrace[,'ss']==1, 'asymm'])
hist(ztrace[ztrace[,'ss']==2, 'asymm'])
dev.off()

pdf(file="wdx.pdf")
hist(ztrace[ztrace[,'ss']==1, 'wdx'])
hist(ztrace[ztrace[,'ss']==2, 'wdx'])
dev.off()

met<-"dNx"
pdf(file=paste(met,".pdf", sep=''))
hist(ztrace[ztrace[,'ss']==1, met])
hist(ztrace[ztrace[,'ss']==2, met])
dev.off()

pdf('corrs.pdf')
pairs(ztrace[,c(1, 3:5)])
dev.off()